Sudden Death of Entanglement induced by Polarization Mode Dispersion 
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We study the decoherence of polarization-entangled photon pairs subject to the effects of polariza- 
tion mode dispersion, the chief polarization decoherence mechanism in optical fibers. We show that 
fiber propagation reveals an intriguing interplay between the concepts of entanglement sudden death, 
decoherence-free sub-spaces and non-locality. We define the boundaries in which entanglement-based 
quantum communications protocols relying on fiber propagation can be applied. 
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Entanglement between particles is a fundamental fea- 
ture of quantum physics. Just as fundamental is the phe- 
nomenon of decoherence that takes place when the en- 
tangled quantum system interacts with the environment. 
One of the most intriguing recent discoveries related to 
decoherence is the phenomenon of entanglement sudden 
death (ESD) QQ. It manifests itself in an abrupt disap- 
pearance of entanglement once the interaction with the 
environment reaches a certain threshold. Beyond 
the interest that it attracts as a fundamental physical 
phenomenon, decoherence plays a central role in quan- 
tum communications. The security of recent quantum 
key distribution protocols explicitly relies on the nonlo- 
cal properties of entanglement, quantified in terms of the 
violation of a Bell-type inequality . Therefore estab- 
lishing the relation between the violation of non-locality 
and ESD, which is very interesting from a standpoint of 
basic physics, has a potentially large impact on the new 
area of quantum communications. 

A configuration that provides an excellent platform 
for the controlled study of decoherence is that of polar- 
ization entangled photon-pairs, distributed over optical 
fibers. In this scheme the main source of decoherence 
is the residual optical birefringence randomly accumu- 
lating along the fiber. While an alternative entangle- 
ment scheme, insensitive to birefringence, has been pro- 
posed [H, [13] , the ease with which light polarization can 
be manipulated using standard instrumentation leaves 
polarization-entanglement the configuration of choice in 
many situations ll|. Moreover, numerous sources of po- 
larization entangled photons suitable for use with stan- 
dard fibers have recently become available Hence, 
understanding the relation between non-locality and ESD 
as well as the ultimate limits imposed by fiber birefrin- 
gence on the distribution of polarization-entangled pho- 
tons in fibers is a problem of utmost importance. 

The fact that optical birefringence is a major polar- 
ization decoherence mechanism has been known for a 
while. Indeed bircfringent crystals have been used exten- 
sively for the creation and manipulation of special quan- 
tum states, such as the MEMS [H [3] 



or Werner states 



13L 15 1 . Similarly, birefringent crystals have also been 
used for the controlled demonstration of decoherence-free 
subspaces H 13] • Yet, the arbitrary birefringence char- 
acterizing fiber-optic transmission, produces a previously 
unobserved combination of physical effects. 

The accumulation of randomly varying birefringence in 
fibers leads to a phenomenon known by the name of po- 
larization mode dispersion (PMD) [lSj. Since the anal- 
ysis of the general case of PMD is quite cumbersome, 
we limit ourselves to the simplest regime of operation in 
which the optical bandwidth of the photons is small in 
comparison with the bandwidth over which PMD decor- 
relates (ljj . In this regime, without loosing the essence 
of the problem, the overall effect of PMD resembles that 
of pure birefringence in the sense that it causes an inci- 
dent pulse to split into two orthogonall y p olarized com- 
ponents delayed relative to each other [18|. The polar- 
ization states of these two components are known as the 
principal states of polarization (PSP) and the delay be- 
tween them is called the differential group delay (DGD). 

In contrast to the controlled environment of (l3j-[l7|. 
both the PSP and the DGD of real fibers vary stochas- 
tically in time. Since typical time constants character- 
izing the decorrelation of PMD in optical fibers are as 
long as hours, days and sometimes months (joj . PMD 
evolution can be considered adiabatic in the context of 
quantum communications protocols. Thus the density 
matrix describing the quantum state needs to be evalu- 
ated as a function of the arbitrary values of the instanta- 
neous PMD. As a consequence, the parameters of interest 
obtained from the so evaluated density matrix are also 
PMD dependent. The temporal statistics for those pa- 
rameters could be in principle determined by application 
of the proper PMD statistics. An approach, in which 
the randomness of PMD is accounted for in the density 
matrix itself 2l|, implicitly assumes ultrafast PMD dy- 
namics and leads to fundamentally different results. This 
previously unstudied reality, produces important conse- 
quences to the dynamics of decoherence between polar- 
ization entangled photons. 

In this letter we formulate a quantitative approach to 
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studying PMD-induced disentanglement. We consider 
the evolution of an arbitrary two-photon state maxi- 
mally entangled in polarization as each photon propa- 
gates through a fiber with PMD. Our studies demon- 
strate that the unequal and increasing differential delays 
in both arms always lead to entanglement sudden death 
P, |i| for all but two special PSPs orientations. That is, 
our channel causes an abrupt drop to zero in concurrence 
while a single photon subjected to the same environment 
depolarizes asymptotically. Contrary to that, when both 
delays are sufficiently close in value there exists a range 
of PSP orientations for which concurrence does not van- 
ish. This is related to the existence of decoherence free 
subspaces and offers an opportunity for non-local PMD 
compensation. Finally, when only one photon experi- 
ences PMD, the concurrence decays gradually for every 
PSP orientation. Besides concurrence, we calculate the 
S parameter of the Clauser-Horne-Shimony-Holt Bell's 
inequality. Comparing the loss of entanglement (concur- 
rence C = 0) and violation of locality (S > 2) we discover 
an intriguing empirical relation between them. 

We assume a source in which a pair of polarization 
entangled photons is generated either via spontaneous 



parametric down conversion 22], or by using four-wave 
24] . In either case, the quantum state of the 
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generated photon-pair can be expressed as 
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where the left-hand side of the tensor product represents 
the frequency waveform, whereas the right-hand side rep- 
resents polarization modes. The frequency variable uj 
denotes the offset from the central frequency, which is 
equal to the pump-frequency in sources relying on four- 
wave mixing and to half the pump frequency in sources 
based on spontaneous parametric down conversion. The 
function /(w) represents the effect of phase- matching as 
well as the possible effects of filters, as we shall see be- 
low. Notice that normalization of the state \tp) implies 
that J dw|/(w)| 2 = 2ttT~ 1 with T being the integration 
time of the detectors used in the set-up. In what fol- 
lows, the state (JT|) will be referred to as the input state 
of the system, which consists of the two optical paths 
that lead the entangled photons from the source of en- 
tanglement towards its users, conventionally referred to 
as Alice and Bob. The terms u A B and v/ A B are the 
Jones vectors that correspond to the excited polarization 
states of Alice's and Bob's photons, respectively. We use 
primes to denote orthogonality in polarization space, so 
that u A B ■ u' A B — 0. The phase factor a is introduced 
for consistency with the experimental generation of po- 
larization entangled photon pairs [23|. Notice that the 
frequency dependent part in (JT]) can be re-expressed as 
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with /(r) being the inverse Fourier transform of f(uj). 
As can be deduced from the form of Eq. ([2]), |/(r)| 2 
is proportional to the probability density function that 
Bob's photon precedes Alice's photon (or vice versa) 
by r. For brevity, we will denote the polarization 
dependent part in the tensor product (jTJ) by \ip p ) — 
[\u A ,u B ) + e ta \u' A , u' B )] /V%, whereas the expression in 
Eq. ([2]) will be shortly denoted as \f(t A — t B )). The over- 
all state is then expressed as \tp) = \ f(tA — t B )) <£> \ip P )- 

Let us now represent the state in terms of the prin- 
cipal states of the PMD in the two arms. We denote by 
{s A ,s' A } and {s B ,s' B } the pairs of Jones vectors that cor- 
respond to the PSP along the paths of photons A and B, 
respectively. We now represent \tp p ) in the basis of the 
PSP modes as follows 

\%) = m{k A ,s B )+e^\s' Al s' B ))/V2 

+ mi\s A ,d B )-^W A ,s B ))/V2, (3) 

where the coefficients 771 and r\i are given by 

Vi = (s.a ■ U.a) (s.b ■ Mb) + e iQ (s A ■ u A ) (s B ■ u B ) (4) 
m = (s A ■ u A ) (s' B ■ u B ) + e iQ (s A ■ vf A ) (s' B ■ u! B ) (5) 

and where a,; is defined through the relation 7^ = 
|»7i| exp (i(a — &i)/2). Also, note that, as is implied by 



(1) state normalization, \rji\ 



\m\ 



1. The quantity |ryi| 



is related to the alignment between the input two-photon 
state ([1]) and the PSP. Thus, for example, in the case of 
1L A = s A and u B = s B , the value of |?7i| 2 is unity. In the 
presence of PMD the arrival time of the A-photon is de- 
layed by t a /2 in the s A polarization and advanced by the 
same amount in the s' A polarization, and the B-photon 
undergoes a similar process. Therefore, the output state, 
i.e. the two-photon state after propagating through me- 
dia with PMD, can be expressed as 



dt A dt B f(t A -t B )\t A ,t B ) (2) 



= ^\f(t A -t B - T -^^))®\s A ,s B ) 

+ ^\f(t A -t B -^±^))®\s A ,s! B ) 

- ^\f(t A -t B + ^±^))®\s> A ,s B ) 

+ r ^-\f(t A -t B + T -^^-))®\i A ,d B ). (6) 



The density matrix that characterizes the detected field 
is given by p = J dt' A dt' B (t' A , t' B \ip out ) (ipout\t' A , t' B ), where 
tracing over the time modes is performed in order to 
account for the fact that the photo-detection process 
is not sensitive to the photon's time of arrival (within 
the detector's integration window). The elements of 
the resulting density matrix, establishing the correspon- 
dence (s A ,s B ) <H- 1, (s' Al s B ) <H> 2, (s A ,s' B ) ■(-> 3, and 
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FIG. 1: Concurrence versus ta and j»7i| 2 . In (a) tb = ta, 
whereas in (b) tb — 1.7B~ . Light-green: S > 2. Dark- 
green: S < 2. Red: C = and S < 2. 



(s!a,s' b ) O 4 are then given by 
Pit = P44 = | 2 /2 

P22 = P33 = | r? 2 1 2 /2 

P3i = -P42 = vlmRf ( t b) /2 

P21 = -P43 = -VlV2 e ' aR f ( T A) /2 
P41 = Wfe'^^-TB)^ 

P23 = -(v*2) 2 e* a R f (T A +T B )/2, (7) 

where Rf(r) = T J dtf*(t)f(t + r) is the autocorrela- 
tion function of /(f), normalized such that Rf(0) = 1. 
The function /(w) which defines the frequency con- 
tents of the two generated photons accounts for the 
phase matching spectrum, as well as for filtering ap- 
plied to the two generated photons. In this case f(u>) = 
/pm(w)i? J 4(o;)if j B(— lu), where Ha(uj) and Hb{w) denote 
the transfer functions of Alice's and Bob's filters, respec- 
tively, and where f pm (ui) represents the phase-matching 
spectrum. In most applications, the filters are much 
narrower than the phase matching spectrum, in which 
case f pm (ui) can be replaced by a constant, such that 
f(oj) oc H A (u)H b (-oj). Notice that the effect of PMD 
scales with the width of the autocorrelation function 
Rf(r), which is in turn determined by the overlap band- 
width of Alice's and Bob's filters, namely by the width 
of \Ha{uj)H b {— w)| 2 . For illustrative purposes, in all the 
numerical examples considered in what follows, we will 
assume that \H a {oj)\ 2 and \H B (co)\ 2 are Gaussian func- 
tions of root mean square bandwidth B and central fre- 
quencies oja an d Wb = — lo a respectively. 

We now turn to the characterization of the degree of 
entanglement of the PMD-affected two-photon state. In 
the presence of PMD, tracing over the time of arrival puts 
the system in a partially mixed state. The extent of this 
process can be quantified with the help of several pro- 
posed entanglement metrics 25l427j . We choose to calcu- 
late purity, concurrence and Bell's S parameter, thus 
assessing the largest possible violation of Bell's inequal- 
ity in the CHSH definition 28[. Because the individual 
photon states are maximally mixed, the two-photon den- 



sity matrix can be reduced to a Bell-diagonal form by a 
proper change of basis [2^] . This enables a fully analyt- 
ical evaluation of C and S. Note that a Bell-diagonal 
matrix is defined by three real parameters only. In the 
case of PMD they are ta, t b and |?7i| 2 and the functional 
dependence of C and S on them is given in [2^] . 

In the simplest case of PMD present in only one of 
the two fibers, as described, for example, by t b = 0, 
the concurrence is given by C = \Rf(T A )\; in this case 
the concurrence is independent of the PSP orientation 
and can only decay asymptotically with ta- Correspond- 
ingly, the S parameter acquires the maximum value com- 
patible with such concurrence, that is S = 1\J\ + C 2 
30], indicating unconditional violation of Bell's inequal- 
ity. Note that the two cases |77i| = 1 and = 0, 
where concurrence simplifies to C — \Rj(ta — tb)\ and 
C = \Rf(jA + t b )\ respectively, are equivalent to that 
of single- arm PMD, with corresponding nonzero DGDs 
equal to ta — t b and ta + t b . Remarkably, for ta = t b 
and 1 771 1 = 1, the concurrence is unity, regardless of the 
DGD magnitude. This result is quite interesting and it 
is directly related to the concept of decoherence-free sub- 
spaces 2, 12 1- In this situation the output state, when 
expressed in the basis of the PSP, is a superposition of a 
state in which both photons are delayed, with a state in 
which they are both advanced (see Eq. ([3])), such that 
they reach the detectors simultaneously. Since in this 
state, knowledge of the photon's times of arrival discloses 
no information on their polarization states, tracing out 
time involves no loss of information. Decoherence-free 
subspaces would not be allowed if PMD dynamics were 
fast on the scale of measurements, as assumed in [21 [. 

The dependence of the two-photon state decoherence 
on the PMD parameters in the general case is more cum- 
bersome, as it is governed by the two DGD values ta and 
tb and by the PSP orientation, accounted for by \r}\ \. For 
illustrative purposes, we plot in Fig. [T]the concurrence 
as a function of ta (normalized to £? _1 ) and \r]i\ 2 for two 
different settings of tb , so as to describe PMD effects for 
the most relevant realizations of PMD parameters. In 
Fig. QJi we plot the concurrence for the case of identical 
DGDs, t b = ta, whereas in Fig. Q}) the value of t b is 
fixed and equal to 1.7B -1 . The range of values for which 
entanglement disappears entirely (i.e. C = 0) is empha- 
sized in red color, whereas the green colored regions of 
the surface correspond to settings in which entanglement 
exists (i.e. C > 0). In the light green area S > 2, whereas 
in the dark green area S < 2. In the case where the DGD 
values are equal (Fig. la), the concurrence approaches 
unity when \i]i\ 2 — > 1. That is because in this situation 
the input state is given by the first term in Eq. ([3]), and it 
is not affected by the loss of time of arrival information. 
In the opposite limit, when |rii| = 0, the concurrence is 
given by C — \Rf(2r A ) \ and reduces towards zero asymp- 
totically for large DGD values. Figure [Tp also illustrates 
that concurrence can be unity only when the DGD values 
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FIG. 2: ESD probability (see text) versus ta and tb normal- 
ized to The grey line is the boundary S = 2. Below 
and to the left of the grey line S > 2 for all values of |?7i| 2 , 
whereas to its right S < 2 for some range of \rji\ 2 . 



in the two arms are equal. 

The most interesting feature in Fig. [T] is the abrupt 
transition of the concurrence to zero when either ta, 
or |?7i| 2 are varied continuously. The abrupt decay of 
concurrence is contrasted with the asymptotic decoher- 
ence that would be experienced by a single photon if it 
evolved in the same environment. Indeed, the purity of 
a polarized single-photon pulse characterized by a Jones 
vector u and a time-mode distribution g(t), transmit- 
ted through a fiber with DGD r and PSP s, is given by 
P = 1 - % ' s\ 2 (l-\u- 



the autocorrelation function of g(t). This is a manifes- 
tation of the entanglement sudden-death (ESD) that has 
been previously reported in other physical systems [lj-Q . 

In general, PSP orientation in optical fibers varies 
faster than the DGD, resulting in a uniform distribution 
of the parameter \t]i\ 2 between and 1. For each given 
combination of ta and tb, we evaluate the fraction of 
the interval between and 1 in which ESD occurs. This 
quantity, which can be interpreted as the probability of 
ESD (conditioned on the DGD values) is illustrated in 
Fig. 2. Note the existence of a completely white re- 
gion that shows the range of DGD values in which en- 
tanglement does not disappear for any value of 1 771 1 2 . In 
contrast to that, the darker tones mark areas for which 
ESD occurs for some range of \r]i\ 2 . The color progres- 
sively turns dark for highly differing DGD values. On the 
other hand, when increasing DGDs remain nearly equal, 
the probability of ESD is about 0.5 

Many applications involving entanglement rely directly 
on the violation of Bell's inequality 0, [1] and, there- 
fore, we compute the maximum value of the CHSH S- 
parameter [28[ for the density matrix Eq. ([7]) 31]. In Fig. 
1, the range of parameters in which S > 2 — meaning 
that the CHSH inequality is violated — is colored by light 
green, whereas the dark green part of the surfaces rep- 
resents states in which C > 0, but there is no violation 
of the CHSH inequality (S < 2). The thick gray line in 



Fig. 2 also marks the S = 2 boundary. Below and to the 
left of this line S > 2 for all |^i| 2 values, but to its right 
S may be smaller than 2 for some relative PSP orien- 
tations. Intriguingly, this non-locality boundary (thick 
gray line) nearly perfectly reproduces the boundary to 
the ESD-free region (the white area) if the scale on both 
axes is stretched by a factor of 1.5. 

To conclude, we have carried out what we believe to be 
the first quantitative analysis of the decoherence of polar- 
ization entangled photons propagating in optical fibers. 
Our study shows that PMD leads to entanglement sud- 
den death in all but a well defined restricted set of re- 
alizations of fiber birefringence. In addition, we demon- 
strated how decoherence-free subspaces can be reached 
via non-local PMD compensation. The ultimate limits 
imposed by fiber birefringence to applications based on 
non-local properties of polarization entanglement were 
shown to be intriguingly related with the phenomenon of 
entanglement sudden death. 
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In this supplement we derive analytical expressions for the purity, concurrence and S pa- 
t^- ' rameter of the Clauser-Horne-Shimony-Holt-type Bell's inequality, which characterize a pair 

of polarization entangled photons transmitted through two optical fibers with polarization 
mode dispersion. 
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d 1 1. INTRODUCTION 

P 

cr . 

In our Letter [l( we derive the density matrix describing the polarization state of two 
polarization-entangled photons transmitted over optical fibers with polarization mode dispersion 
(PMD). The density matrix of Eq. (7) in [l| (we will further down refer to this equation as Eq. 
(M7)) is expressed in the frame of reference of the principal states of polarization (PSP), in which 
the derivation is greatly simplified and the final result appears quite compact. However, the an- 
alytical evaluation of some of the parameters characterizing the entanglement of the two-photon 
state is more convenient in a different frame of reference, in which the density matrix appears in 
the so-called X-shape [1, 0] ■ 

In the first section of this supplement we study the purity of the two-photon state. In the 
following sections we provide analytical formulae for concurrence and CHSH 5-parameter. 

• i— i 

X 

H 2. PURITY OF THE TWO-PHOTON STATE 

The purity V of the two-photon state can be evaluated through straightforward algebra from 
Eq. (M7) 

V = |7?l|4 + |7?2|4 +ki| 2 M 2 [|^(r A )| 2 + |^(r B )| 2 ] 

+ l[\ m \±\R f (T A -T B )\* + \ m \*\R f (T A +T B )\ 2 ]. (1) 

Obviously, the purity is equal to 1 when ta = tb = 0, that is when there is no PMD at all. 
Yet interestingly, it remains equal to 1 when \rji\ 2 = 1 and the DGD values are equal, re- 
gardless of their magnitude. This property is a manifestation of the decoherence-free subspace 
concept. On the other hand, as the DGD values become larger than the inverse overlap band- 
width of Alice's and Bob's filters B^ 1 (see [jij for the definition of B), the purity approaches 
|t7i| 4 [l + |.R/(At)| 2 ] /2 + 1 772 1 4 /2 , with At = ta — t~b- Then, the minimum purity with respect to 
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r/i and rj 2 is [l + \Rf(Ar)\ 2 ] / [2 + ^/(Ar)] 2 ] /2 and it ranges between 1/4, for At = 0, and 1/3, 
forAT»r 1 . 



3. EXPRESSION FOR CONCURRENCE AND 5- PARAMETER 

In this section we derive the expressions for the concurrence and the S parameter of the CHSH 
version of Bell's inequality in terms of the parameters characterizing the PMD of the two fiber 
arms. This goal is most conveniently achieved by first expressing the density matrix in a basis 
where only the terms on its two diagonals are not zero, a form known as X-shape representation 
@, 0] . Subsequently the procedure requires defining the matrix S whose components are given by 
Sum = Trace (po n (g) cr m ), where a n (n = 1, 2, 3) are the Pauli matrices. Then the concurrence and 
the S parameter are expressed in terms of the eigenvalues of S T S, where T denotes transposition. 
Eventually, the eigenvalues of S T S are related to the PMD parameters. 



3.1 The X-shape representation of the two-photon density matrix 



Since the individual photons are completely depolarized, the two-photon state belongs to the 
class of the Bell-diagonal states fl and hence to that of the so-called X-states 0, 3 . Namely, there 
exists a representation in which the density matrix takes the form of 



pn p u 
P22 P23 

P32 P33 

Pai p 44 

where all entries are real with pu = p^A, P22 = P33 = 1/2 - 
representation concurrence is given by the formula [3]: 



(2) 



Pn, Pu = Pai and ,932 = p 2 3- In this 



C = 2max{0, |p 23 | - pu, \p u \ - P22] 



(3) 



The representation in which the density matrix is in the form of Eq. ([2]) is constructed as follows. 
Using the coefficients s nm = Trace (pa n <g> a m ) which form a 3 x 3 real matrix S, we find the 
eigenvectors of S T S and of SS T , which could be viewed as vectors representing states of polarization 
in the Stokes space 0] • One then chooses the two orthogonal Jones vectors that correspond to one 
of the eigenvectors of S T S to span the subspace corresponding to the first photon. Similarly an 
orthogonal pair of Jones vectors corresponding to one of the eigenvectors of SS T is chosen to span 
the subspace corresponding to the second photon. The overall space, which is the tensor product 
of the two subspaces, is then spanned by the two pairs of Jones vectors that were selected 0]. In 
this process, the eigenvector of S r S and the eigenvector of SS r that the desired representation 
is based upon, are selected arbitrarily. Moreover, it is irrelevant which of the two pairs of Jones 
vectors is assumed to represent a specific photon's subspace, namely the two pairs can always be 
interchanged without compromising the process. 

In this representation one can verify, using Eq. ([2]) along with the definition of the s nm elements, 
that S is diagonal and therefore 



S T S 



(pu + P23) 2 







(pu ~ P23) 









(p 



11 



P22) 



(4) 



This form will be later used for relating the concurrence and the S parameter to the quantities 
characterizing PMD through the eigenvalues of S T S. 
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3.2 Expressing concurrence in terms of the eigenvalues of S T S 

The aim of this section is to demonstrate the following expression for concurrence in terms of 
the eigenvalues of S T S: 

C= imax{o,Ai /2 + A2 /2 + A3 /2 -l}, (5) 

where we denote these eigenvalues as Ai, A2 and A3. For the calculation that follows, it is convenient 
to assume that the eigenvalues are sorted such that Ai > A2 > A3. 

Our first goal is to express the Bell-diagonal density matrix in terms of Ai, A2 and A3. Let us 
consider a basis for subspace A such that the corresponding Stokes unit vector b& is an eigenvector 
of S T S belonging to the largest eigenvalue, namely S t S6a = XibA- Consider also a basis for 
subspace B such that the corresponding Stokes unit vector bs is an eigenvector of SS T defined 
though bs = S6j4./| I Sbv4 1 1 , where || • || denotes vector length. 

We denote the correlation between the measurements of coincident photons along the po- 
larizations defined by the Stokes vectors 6^ and &b as E(bA,bB) = P c — (I — P c ), where 
Pc = P11+P44 = 2pn is the probability of a coincident detection of photons A and B passing through 

polarizer beam splitters aligned to %a and bs- From Eq. (7) of jij, E(bA,bs) = • S&b = A^ 2 , it 

1/2 

follows that p\\ = P44 = (1 + A x )/4. Also from the condition of unit trace of the density matrix, 

P22 = Pzz = (1 — A^ 2 )/4. This yields Ai as the third diagonal element of S T S. 

The remaining off-diagonal terms of the density matrix follow from equating the first and the 
second diagonal elements of S r S to A2 and A3 (in some arbitrarily assumed order). This yields two 

possible solutions: either |p23 1 = \PZ2\ = A^+Ai!/ 2 and \p\^\ = \pa\ = A<|/ 2 — Aj!/ 2 , or \p2%\ = |/o>32 1 = 

1/2 1/2 1/2 1/2 

A 2 — A 3 and \p\4\ = \pa\ = A 2 + A 3 . However, only one of them is compatible with nonzero 

1 /2 1 /2 

values of concurrence. In fact, using Eq. ([3|) for the choice \p2z\ = \PZ2\ = A 2 + A 3 and \p\4\ = 

I, vl/2 ,1/2 . , , n f n x l/2 . >l/2 .1/2 1 ,1/2 ,1/2 ,1/2 A ic, n 

I/O41I = A 2 ' —A 3 yields G = max<0, — \{ + A 2 + A 3 ' — 1, X{ + A 2 — A 3 ' — 1 > /2 = 0. 
Indeed, since all eigenvalues are less than unity and Ai is the largest, the following inequality is true, 
— A^ + A^ + Aj!/ 2 — 1 < 0. On the other hand, the positivity of p implies that A 2 ^ 2 — Ag^ 2 < 1 — Aj^ 2 , 
which yields Aj^ 2 + A^ 2 — A^ 2 — 1 < 0. Hence the correct solution is |>c?23 1 = \pm\ = A^ 2 — A^ 2 
and I/014 1 = \pu\ = A 2 ^ 2 + A^ 2 and the density matrix reads as: 



P= 4 



1 + A} /2 C(A2 /2 + Ag /2 ) 

1-Ai /2 e(A2 /2 -Ag /2 ) 

^Ap-Ag 72 ) 1-Ai /2 

C(A 2 /2 + A3 /2 ) 1 + Ai /2 



(6) 



where C and £ can arbitrarily be set to ±1. This indetermination reflects the phase ambiguity in 
the conversion from Stokes to Jones representation of a polarization state 0]. 

Finally, applying the simplified concurrence formula Eq. ([3]) to the density matrix in Eq. ([6]) 
yields the expression for concurrence given in Eq. ©. 



3.3. Expressions for the eigenvalues of S T S in terms of the PMD parameters 

The analytical derivation of the eigenvalues Aj requires to solve a third order algebraic equation. 
This task is in general quite cumbersome. However, for the density matrix in the shape of Eq. 
(M7) this is possible under two conditions that make the matrix real. First, we make a realistic 
assumption on the symmetry of Alice's and Bob's filters, \Ha(uj)Hb{—uj)\ = \Ha(—w)Hb(uj)\. 
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Under this assumption Rf(r) is real. Secondly, we perform a change of basis. In particular, 
we use a basis (\s A ),\s' A )) for the subspace A where \s A ) = e lipA / 2 \s A ) and |s^) = e~ lifA ^ 2 \s' A ) 
with ip A = a + Im[log(77i7/2)]. Similarly, we use a basis (\s B ) , \s' B )) for the subspace B where 
\§_b) = ^ b ^\§.b) an d \s'b) = e ~ lipB ^ 2 \dB) w hh ifB = Im [log (771 7/2)]- In this representation the 
density matrix of Eq. (M7) can be simplified by setting a = and by replacing 771 and 772 with 
their absolute values |7?i| and 1 772 1 , thus making the matrix real. In this representation the matrix 
S T S appears as: 



S T S 



mn 777-13 
777,22 

mi3 77733 



(7) 



with 



77711 



[\m\ 2 Rf(TA-T B )-\v2\ 2 Rf(TA + T B )] +4\ m \ 2 \ m \ 2 R}{T B ) (8) 

"122 = [|77l| 2 i?/(TA - TB) + \rj 2 \ 2 Rf(TA + T B )} 2 (9) 

77733 = 4|7 ? i| 2 |7 ?2 | 2 J Rf(r j4 ) + (2|7 ? i| 2 -l) 2 (10) 

77713 = -2[\m\ 2 Rf( T A- t b ) - \m\ 2 Rf(TA + T B )]\m\\m\Rf(TA) 

+ 2\ri lV2 \R f (T B )(2\ m \ 2 - 1) (11) 
The eigenvalues of S T S can be evaluated with the following result 

il = (mn + 777 33 )/2 + [(77711 + 777 33 ) 2 /4 - (7771177733 - 777 2 3 )] V2 (12) 

£ 2 = 77722 (13) 

1 /2 

h = (mn + m 33 )/2 - [(mn + m 33 ) 2 /4 - (mnm 33 - m^ 3 )] , (14) 
where £i, £2 and £3 are the unsorted eigenvalues of S T S. The concurrence is then given by 



c 2 



im^{0,<! /2 +^ /2 +4 /2 -l}. (15) 



The largest violation of CHSH inequality is given by S = 2 V / Ai + A2, where Ai and A2 are the two 
largest eigenvalues of S r S [6]. From Eqs. (12) and (14) it follows that l\ > £3, and therefore 

S = 2v^i +max{e 2 ,e 3 }. (16) 

Equations (115p and fjl6[) are the main result of this supplement. They give the expressions for 
largest violation of Bell's inequality and concurrence as a function of the instantaneous PMD 
values. 



3.4 Entanglement sudden death threshold for large differential group delays 

We conclude this supplement by evaluating the concurrence of the output two-photon state in 
a particular case of interest, that is when the two DGDs are much larger than the inverse photon 
bandwidth. In this regime one can derive the following expression 

C = max {0, [1 + \R f (AT)\] |t7i| 2 - l}, (17) 

where At is the difference between the two DGDs. This results can be used to characterize a 
threshold seen in Fig. la of the Letter. In that figure it appears that for two equal DGDs there is a 
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threshold value of \rji\ 2 = 0.5, above which ESD does not occur for increasing DGDs. In fact, such 
a threshold exists as long as the difference between two DGDs is finite. Using the above asymptotic 
expression for concurrence this threshold can be estimated as l^i^hl 2 = [1 + |-R/(At)|] 
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